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ABSTRACT 



0^ Aims. We investigated the physical properties and dynamical evolution of Near Earth Asteroid (NEA) (190491) 2000 



FJio in order to assess the suitability of this accessible NEA as a space mission target. 

Methods. Photometry and colour determination were carried out with the 1.54 m Kuiper Telescope (Mt Bigelow, USA) 
and the 10 m Southern African Large Telescope (SALT; Sutherland, South Africa) during the object's recent favourable 
apparition in 2011-12. During the earlier 2008 apparition, a spectrum of the object in the 6000-9000 Angstrom region 
was obtained with the 4.2 m William Herschel Telescope (WHT; Canary lis, Spain). Interpretation of the observational 
results was aided by numerical simulations of 1000 dynamical clones of 2000 FJio up to 10 6 yr in the past and in the 
.£j ■ future. 

Results. The asteroid's spectrum and colours determined by our observations suggest a taxonomic classification within 
the S-complex although other classifications (V, D, E, M, P) cannot be ruled out. On this evidence, it is unlikely to 
be a primitive, relatively unaltered remnant from the early history of the solar system and thus a low priority target 
for robotic sample return. Our photometry placed a lower bound of 2 hrs to the asteroid's rotation period. Its absolute 
magnitude was estimated to be 21.54±0.1 which, for a typical S-complex albedo, translates into a diameter of 130±20 
1 m. Our dynamical simulations show that it has likely been an Amor for the past 10 5 yr. Although currently not Earth- 

qq _ crossing, it will likely become so during the period 50 — 100 kyr in the future. It may have arrived from the inner 

or central Main Belt > 1 Myr ago as a former member of a low-inclination S-class asteroid family. Its relatively slow 
■ rotation and large size make it a suitable destination for a human mission. We show that ballistic Earth-190491-Earth 

transfer trajectories with AV < 2 km s _1 at the asteroid exist between 2052 and 2061. 

Key words. Minor planets, asteroids: individual: 2000 FJio - Methods: observational - Methods: numerical 

1. Introduction size, shape, surface roughness, rotational state and spectral 

type - ranges from poor to non-existent. Part of the prob- 
The population of near-Earth asteroids (NEAs) contains a lem stems from their Earth-like orbits; slow keplerian shear 
; small fraction of objects in low inclination, low-eccentricity generally places them beyond the reach of Earth-based ob- 
. orbits similar to the Earth's. These NEAs are considered ser vatories except during the few months that they spend 
attractive targets for in situ investigation by robots or hu- in proximity to the Earth every decade or so. The situation 
mans ([Hasegawa et al|2008|;|Abell et al|2009|;|Mlchel et al.| is also not helped by their small sizeSj typically a few tens 
120PJ; |Lauretta et al.||2qi0t |Elvis et al.||201l|). However, the of metr es, so even when near the Earth their study is the 
attractiveness of individual objects as targets for either exclusive purvue of large-aperture instruments. Primitive 
robotic or human missions is mired by the currently poor NEA taxonomies - those belonging to classes B, C, D, and 
knowledge of their orbits and physical properties. Many p _ are pre ferred as mission targets for robotic sample re- 
have been observed only on a single apparition resulting turn as they are thought t o be relatively unal tered relics of 
in large projected uncertainties in their future position. t he early solar system (ee lMichel et"alf [2009). In addition, 
In addition, knowledge of properties that are important fast ro tators are unsuitable as mission targets due to the 
from an operational as well as a scientific standpoint - added complexity of operations in close proximity to such 
objects. 
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This paper reports on a study of NEA (190491) 2000 
FJio, an object that is accessible from the Earth (Section 
[2]). A programme of observations of 190491 was carried out 
with the 4.2 m William Herschel Telescope (WHT), the 10 
m Southern African Large Telescope (SALT) and the 1.5 m 
Kuiper telescope (Section [3]) aiming to constrain its taxo- 
nomic type, size, and rotational state (Section 2|) . Those 
were combined with numerical simulations of 190491's or- 
bital evolution (Section [5]) to provide context for the ob- 
servational characterisation and help us draw conclusions 
on the object's likely origin (Section [5]). Its accessibility 
from the Earth was quantified by constructing direct, two- 
way keplerian trajectories between the Earth and the as- 
teroid (Section [7]) . A summary of our findings is provided 
in Section [51 

2. The Asteroid 

(190491) 2000 FJio was discovered by the Spacewatch sur- 
vey on 25 March 2000. Based on its orbital parameters 
(a = 1.32 AU, e = 0.23, i = 5°) it is classified as 
an Amor Near Earth Asteroid (NEA). Its absolute mag- 
nitude H — 20.9 implies a diameter ranging from 110- 
390 m for an albedo range 0.05-0.5. Its Earth Minimum 
Orbit Intersection Distance (MOID) is 0.055 AU, slightly 
higher than the threshold of 0.05 AU for classifying it as 
a Potentially Hazardous Asteroid (PHA). It is one of the 
most accessible spacecraft targets, ranking 124th out of 
8857 object entries in the list of Near Earth Asteroid AV 
for spacecraft rendezvouses of May 2012. Its AV as stated 
in that list is 4.567 km s _1 , slightly above the threshold 
that would classify it an Ultra Low Delta V (ULDV) ob- 
ject ([Elvis et al.ll201ll ). However, its absolute magnitude is 
the second-brightest within those objects that precede it, 
the first being the ULDV NEO (89136) 2001 US16 with 
H = 20.2 and AV = 4.428 km s _1 . 

3. Observations and Data reduction 

3.1. SALT 

Located at SAAO in South Africa, the 10 m SALT is based 
on the Hobby-Eberly Telescope (Texas, USA) with a pay- 
load moving during th e observations above the stationary 
10 m spherical mirror (jKwiatkowski et aI1l2010L and refer- 
ences therein). It can access objects in the declination range 
from S = —75° to 5 = +10° when they enter the annular 
region on the sky located between zenith distances from 48° 
to 59°. SALT works solely in the queue-scheduling mode, 
in which the exact time of the observations is not known in 
advance. 

Because of its construction, during observations the 
telescope's pupil continuously changes making it impossi- 
ble to perform all-sky photometry. For the same reason, 
twilight flat fields cannot be used in the photometric re- 
duction. Instead, night sky flat fields derived directly from 
the science frames have to be used. 

We observed 2000 FJi with SALT during September 
2011, on the first month of the facility's normal opera- 
tion after an extended period of commissioning. The instru- 
ment of choice was the SALT imaging camera (SALTICAM; 

1 http://echo.jpl.nasa.gov/~lance/delt a-v/| 
delta_v.rendezvous.html 



lO'Donoeue et al J 12003). a mosaic of two CCDs, each with 
two readout amplifiers. 

The aspect data and the observing log are provided in 
Table [TJ During that period the asteroid was located within 
a star field covered by the Sloan Digital Sky Survey (here- 
after SDSS), which allowed us to calibrate our photometry 
with SDSS standard stars. Under normal operations (the 
case for the observations reported here), there is an effort 
to keep target and comparison sources on the same chip 
so as to ease data analysis. This limited our Field of View 
(hereafter FoV) but still left us with enough SDSS stars for 
comparison. 

In order to determine colour indices we used the Sloan 
g, r and i filters. Because of the significant fringing in the 
infrared, the z filter was not used. As the lower bound of 
the range of the possible effective diameters of 2000 FJio 
wa s 110 m, its rotation period cou ld be as short as 10 min 
(e dHergenrother fc Whitelevll201ll ). To check that, the first 
run on 15 Sep, which lasted 20 min, was done with the r 
filter only. A preliminary reduction showed all photometric 
data points to be within the ±0.05 mag range. The next 
run on 22 Sep lasted 50 min. It was executed with all three 
filters, in the following sequence: 10 x r, 3 x g, r, 3 x g, 
r, 3 x g, r, 3 x i, r, 3 x i. Sequencing the observations in 
this way allowed us to use r exposures to search for bright- 
ness variations, and - if necessary - reduce the magnitudes 
obtained in the other two filters to the same reference level. 

Data reduction was done in two stages. First, the CCD 
frames were corrected for cross-talk and bias usi ng the spe- 
cialised PyjiALTjgackageJ^^ Next, fol- 
lowing iKniazev fc Vaisanenl ([20111) . a flat- field correction 
was carried out using night sky flat-fields obtained from 
the science frames themselves. 

An illumination pattern for each frame was created by 
removing the stars with a median filter and fitting a polyno- 
mial to the remaining background sky. The original frames 
were then divided by these flat-fields. During the obser- 
vations the telescope was dithered every 3 exposures so 
that the stars did not occupy the same positions within 
the FoV. This allowed us to use all frames in a given filter 
- already corrected for the low frequency pattern - to pro- 
duce a global, second-stage flat-field by median combining 
the images. The obtained flat-field mapped the constant 
pixel-to-pixel variations within the FoV. Finally, all images 
in a given filter were divided by the global flat-field frame. 

The instrumental magnitudes of the SDSS stars and of 
the asteroid were then measured with standard aperture 
photometry. Both the flat-field corrections and the photom- 
etry were done with the STARLINK package, maintained 
at the Joint Astronomy Centre of Mauna Kea Observatory. 

A problem in calibrating the instrumental magnitudes 
to the SDSS standard scale was the unknown transforma- 
tion between the photometric systems of SALT and SDSS. 
While this would require a detailed analysis, we instead per- 
formed a simple, first-order check with the available data, 
obtained through the SALT g$, ?s, and is filters . It re- 
vealed some linear trends of rs — r vs. r — g, gs — g vs 
r — g, is — i vs g — i, which led to ~ 0.04 mag differ- 
ences between red and blue stars. Because of this fact, for 
determination of asteroid colours we selected only the so- 
lar type SDSS stars (having 0.4 < g — r < 0.5 mag and 
0.1 < r — i < 0.2 mag), for which the discrepancies be- 
tween both SALT and SDSS photometric systems could be 
neglected. The asteroid magnitudes in each filter were mea- 
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sured with respect to 3-4 SDSS solar colour stars, and the 
results were averaged. This way the uncertainties of the 
SDSS catalogue were minimised. 

On 22 Sep, during 47 minutes, the brightness of the 
asteroid, as monitored in our r filter exposures, did not 
show any systematic changes. This allowed us to average 
the measurements, obtained in r, g, and i filters, and use 
them to compute colour indices. The obtained results are: 
g - r = 0.59 ± 0.03 mag, r - i = 0.21 ± 0.03 mag. 

3.2. Kuiper telescope 

The University of Arizona Kuiper 1.54-m telescope located 
near Mount Bigelow in southeastern Arizona was used 
to obtain V- and R-band photometry of 2000 FJio on 9 
separate dates between September 2011 and March 2012 
(Table [lj. On all 9 nights the asteroid was observed in the 
R-band. On the night of 2011 October 19 the asteroid was 
also observed in the V-band in order to measure its V-R 
color. The instrument used was the Mont4K, a Fairchild 
CCD486 4096 x 4097 CCD, with a FoV of 9.'5 x 9.'5 and 
plate scale of 0.28"/pixel when binned 2x2. 

All data were reduced with the IRAF software package. 
The images were bias-subtracted and flat-fielded with twi- 
light and night sky flat images using the CCDRED package. 
The APPHOT package was used to perform aperture pho- 
tometry of the asteroid and photometric standard stars. 
In order to compensate for variable seeing, the average 
FWHM was measured for each image and the photometric 
aperture was set to a radius of 2 xFWHM. Sky background 
was measured with a circular ring aperture of radius 20 pix- 
els and width of 10 pixels. The sky aperture was centered 
on the position of the measured source. The telescope was 
tracked at the rate of motion of the asteroid and the images 
were shifted and co-added on the motion of the asteroid. 

Photometr ic V- and R-band reference stars from 
lLandoltj (|1992|) were observed at multiple airmasses on each 
night in order to determine the photometric zero point and 
extinction coefficient. 

The V- and R-band images were taken in the following 
sequence: 5xR, 5xV, 5xR, 5xV, 5xR and 5xV. A V-R 
color index of 0.48 ± 0.05 mag was derived from the data. 
The cadence allowed for variability due to the rotation of 
the asteroid to be corrected. During the 37 minutes of ob- 
servations, the asteroid steadily increased in brightness by 
~ 0.3 magnitudes. 

3.3. WHT 

Observations were conducted using the Intermediate- 
dispersion Spectrograph and Imaging System (ISIS) 
mounted on the 4.2m William Herschel Telescope, La 
Palma (Table [lj . Light from the optical system was split 
using the 5300 dichroic (blue cut off/red cut on at 5300 
Angstroms) and directed along the red arm of the instru- 
ment. The red arm uses a red-sensitive 4kx2k pixel RED+ 
CCD with anti-reflection coating. The R158R grating with 
a slit width of 1 arcsecond were used producing a dispersion 
of 1.8 Angstroms per pixel. The total usable wavelength 
coverage from the red arm was 5300 to 10000 Angstroms, 
however the S/N degrades rapidly beyond 9000 Angstroms. 

To ensure the NEA remained within the slit, non- 
sidereal tracking was used with the telescope tracking at 
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Fig. 1. Reflectance spectrum of (190491) 2000 FJi ob- 
tained with the WHT. It is similar to the spectra of S type 
asteroids in the SMASS database. 

the apparent rate of motion of the NEO. The position of 
the NEO was monitored throughout each exposure with 
manual corrections to the pointing position applied when 
necessary. In total 1 x 600s + 3 x 900s exposures were 
obtained. 

Image reduction was performed in the usual manner. 
Bias images were combined and subtracted from the sci- 
ence images. Flat fielding was conducted after removal of 
the spectral profile of the tungsten lamp used to obtain 
the flat field images. The spectra we re extracted using the 
IRAF task apall. Optimal extraction (|Horndll986l) was used 
to improve the S/N of the extracted spectrum. Wavelength 
calibration was achieved using CuAr-CuNe arc lamp spec- 
tra observed at the same pointing position to account for 
the effects of flexure. The solar analog 16 Cyg B was ob- 
served to enable the removal of the solar spectrum from 
the NEA spectrum. Atmospheric correction was perfo rmed 
using tabulated extinction functions for La Palma (King 
1985). These curves were scaled to the airmass of the ob- 
servations of the NEA and solar analog and the corrections 
applied to each spectrum. Asteroid (1) Ceres was also ob- 
served and the spectrum extracted in this ma nner was con- 
sisten t with previously published spectra (e.g lBus fc Binzell 
I2002T ) . Spectra were normalized and the solar spectrum re- 
moved. The result for 2000 FJio is shown in Fig. [TJ 

Chi-square d fitting to the Bus-DeMeo taxonomy 

(|DeMeo et al.ll2009h was conducted. The spectra were re- 
sampled to produce reflectances at the wavelengths used to 
define the various taxonomic types. The best fit was to an 
Sq-type, but we note that the 1 micron absorption band has 
not been fully sampled, making a formal definition difficult. 

4. Physical characteristics 

4.1. Rotation period 

As already explained, the specific construction of SALT lim- 
its its typical continuous observing run to 1 hour. Because 
of this, it is difficult to measure brightness variations with 
periods of several hours or longer. In the case of 2000 FJio 
there was some possibility that - due to its small size - 
it is a fast rotator, with a rotation period <C 2 hr. Fig. [5] 
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Table 1. Aspect data and observing log 



Date 


Obs. time 


r 


A 




A 





V 


Mov 


Exp 


Filter 


Telescope 


YYYY-MM-DD 


(UTC) 


[AU] 


[AU] 


n 


[°] 


n 


[mag] 


["/min] 


M 






9008-1 0-08 


nA-^n - f]5-47 

\J L ±.0\J yjo.^ii 


1.072 


0.081 


26.2 


42.1 


-8.8 


16.7 


4.3 


1 x 600 + 

3 x 900 


pec 


WHT 


2011-09-15 


19:16 - 19:35 


1.1674 


0.1799 


24.0 


327.2 


12.5 


18.7 


1.5 


60 


r 


SALT 


2011-09-22 


21:23 - 22:13 


1.1432 


0.1680 


31.4 


324.2 


9.6 


18.7 


1.5 


60 




SALT 


2011-09-26 


03:56 - 04:08 


1.1325 


0.1636 


35.0 


323.1 


8.0 


18.7 


1.5 


60 


R 


Kuiper 


2011-10-19 


03:40 - 04:17 


1.0654 


0.1441 


57.9 


320.7 


-5.7 


18.9 


1.6 


60 


V, R 


Kuiper 


2011-10-21 


03:35 - 03:46 


1.0605 


0.1429 


59.5 


321.0 


-7.1 


19.0 


1.7 


60 


R 


Kuiper 


2012-01-19 


06:34 - 06:44 


1.0871 


0.1647 


47.7 


74.1 


-36.2 


19.0 


3.0 


60 


R 


Kuiper 


2012-01-26 


07:36 - 07:43 


1.1078 


0.1841 


44.3 


82.1 


-31.4 


19.2 


2.8 


30 


R 


Kuiper 


2012-01-27 


06:49 - 06:56 


1.1107 


0.1871 


44.0 


83.1 


-30.7 


19.2 


2.8 


60 


R 


Kuiper 


2012-02-23 


05:16 - 05:27 


1.2016 


0.3032 


40.1 


104.2 


-15.4 


20.3 


1.9 


60 


R 


Kuiper 


2012-02-24 


05:00 - 05:12 


1.2051 


0.3087 


40.2 


104.8 


-15.0 


20.4 


1.9 


60 


R 


Kuiper 


2012-03-28 


04:31 - 04:45 


1.3233 


0.5348 


42.7 


124.1 


-5.2 


21.8 


1.5 


60 


R 


Kuiper 



Note: r and A are the distances of the asteroid from the Sun and the Earth, respectively, a is the solar phase angle, while A and /3 
are the geocentric, ecliptic (J2000) longitude and latitude. In the next column an average brightness V of the asteroid, as predicted 
by the Horizons ephemeris service, is given. Starting from the ninth column, the table gives the asteroid movement on the sky 
(Mov), the exposure time (Exp), and the filters used (the abbrevation Spec denotes spectroscopic observations) 



presents one of the lightcurves obtained during our observa- 
tions. It shows relative brightness variations of the asteroid 
in the Sloan r filter, and of one of the comparison SDSS 
stars. There are no traces of periodicity in these data with 
a peak-to-peak amplitude greater than 0.05 mag. The ob- 
served scatter is probably caused by imperfect flat fielding 
rather than statistical noise. Assuming that a typical aster- 
oid displays a bimodal brightness variation we can conclude 
the rotation period of 2000 FJio P is longer than twice the 
time span covered by our data (P > 36 min). The data 
obtained on another night with the r filter were analysed in 
the same way. As it also showed no discernible light varia- 
tion, we can raise the lower limit for P to 94 min. Of course, 
there is a small probability that at the time of our observa- 
tions the asteroid was visible close to the pole-on view and 
its brightness variations were difficult to detect - despite 
its short period. 

Similar observations performed on 19 Oct with the 
Kuiper telescope revealed a systematic increase of the aster- 
oid's brightness by ~ 0.3 mag during 37 min. This does not 
contradict our SALT observations as from 22 Sep to 19 Oct 
the solar phase angle almost doubled from 31° to 58°, while 
the observer-centred ecliptic latitude changed by 15°. This 
could have led to an increase of the lightcurve amplitude. It 
is also possible that our SALT observations were performed 
during maximum brightness, when the lightcurve for some 
time could be almost flat. 

The fact that on 19 Oct during 37 min the asteroid 
lightcurve did not show any extremum can be used to raise 
further the lower limit for its rotation period from 94 min to 
about 2 hours (under the assumption of a typical bimodal 
lightcurve). We thus conclude that 2000 FJio is not a fast 
rotator. 

4.2. Taxonomy 

From the photometric standpoint, an asteroid's taxonomy 
is usually determined from its colour indices in the Johnson 
UBVRI system. Instead of transforming our Sloan g-r and 
r-i colours to their B-V and V-R counterparts we used an 
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• Asteroid 
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19:17:46 19:20:38 19:23:31 19:26:24 19:29:17 19:32:10 19:35:02 
UTC 

Fig. 2. Lightcurve of (190491) 2000 FJi obtained in the 
Sloan r filter. There are no apparent periodic variations 
during 20 minutes with a peak-to-peak amplitude greater 
than 0.05 mag. The slope in the linear fit is too small to be 
significant. For comparison, brightness variations of one of 
the SDSS stars are also presented. The observed scatter of 
points is caused by the residual effects of the flat fielding 
procedure. 



analogous classifi cation based on the SDSS magnitudes. 
Ilvezic et all (|200lh used 316 spectra obtained in the SMASS 
survey and convolved them with the SDSS response func- 
tions. As a result they were able to identify different taxo- 
nomic classes on the g-i vs r-i domain (Fig. [3]). The position 
of 2000 FJio °n this graph shows it is most probably an S 
type, but V, D, or E, M, P types cannot be ruled out. This 
conclusion is supported by the V — R = 0.48 ± 0.05 mag 
measured with the Kuiper t elescope, which is consi stent 
with an S-type classification (|Tholen fc Barucclll989|) and 
in agreement with the WHT spectrum. 
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Fig. 3. Colour-colour plot for the 316 asteroids whose spec- 
tra w ere obtained by the SMASS Survey, based on Fig. 10 
from llvezic et al.l (|200lf ). Different taxonomic classes are 
presented by different symbols. The colours of 2000 FJi 
(filled circle with error bars) are most compatible with an 
S type classification, though V, D or E, M, P types cannot 
be ruled out. 
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Fig. 4. Least-square fit of the H-G relation to the asteroid 
phase curve. The scatter of the points is partially caused 
by the changing aspect and the lightcurve amplitude. 



4.3. Diameter 

The observations of 2000 FJio over a wide range of phase 
angles allowed us to plot its phase curve, and to derive the 
absolute magnitude as well as the effective diameter. The 
R magnitudes measured with the Kuiper telescope have 
been transformed to V magnitudes using its V-R colour 
index. The r magnitudes measured with SALT have been 
transformed to V magnitudes with the equations given by 
iJester et al.l (|2005l) (for that we also used the asteroid g — r 
colour index). 

Figure H] presents a plot of V(l,l,a) magnitudes ver- 
sus solar phase angle. The apparent magnitude of the as- 
teroid at each observation has been corrected to unit dis- 
tance from the Sun and the E arth. A least squ are fit 
with the standard H-G relation (|Bowell et al.l [l989h gives 
H = 21.54 ± 0.10 mag and G = 0.25 ± 0.08. The quoted 
uncertainties were estimated by assuming error bars in 



Table 2. (190491) 2000 FJm clone classification statistics 
at different stages during the forward and backward inte- 
grations. 



FORWARD 

time from integration start (yr) 



Outcomes 


+10 4 


+10 5 


+5 x 10 5 +10 6 


q>1.3 AU 








17 26 


Amors 


998 


283 


390 376 


Ap olios 


1 


705 


486 425 


Atens 





11 


106 172 


BACKWARD 




time from integration start (yr) 


Outcomes 


-10 4 


-10 5 


-5 x 10 5 -10 6 


q>1.3 AU 





1 


5 18 


Amors 


996 


979 


667 569 


Ap olios 


3 


19 


296 333 


Atens 








31 79 



V(l,l,a) of ±0.05 mag. It is worth noting the large dis- 
crepancy between H = 20.9 given by the MPC (based en- 
tirely on inaccurate magnitude estimates), and our result. 
It confirms a general rule that the MPC absolute magni- 
tudes derived for NEAs are often underestimated by up to 
0.5 mag. 

The obtained value of G = 0.25 is consistent with 
the a verage value G = 0.23 obtain ed for S type aster- 
oids (jLagerkvist fc Maenussonl Il990l) , further supporting 
our taxonomic classification of 2000 FJio as & n S type ob- 
ject. 

Recently, improved albedo estima tes for NEAs in differ - 
cnt taxonomic classes were derived (| Thomas et al J 1201 If ). 
The average geometric albedo for the S complex was 
found to be pv = 0.26^^03 • Using this value with 
our H magnitude , and the classic formula provided by 
iFowler fc Chillemil (|1992ft . we obtain for 2000 FJi an ef- 
fective diameter of D e g = 0.13 ± 0.02 km. 



5. Dynamical Modelling 

5.1. Method 

To characterise the asteroid's recent orbital evolution and 
constrain its origin, we have generated 1000 dynamical 
clones of 190491 by applying the formal state covariance 
of that asteroid for Julian Date 2455600.5 downloaded 
from As tDy^ to a six-dim ensional gaussian random vec- 
tor (see iDuddv et al.l 12012. for details). The clones were 
then integrated 10 6 yr in the past and in the future un- 
der an 8-planet model of the solar system. The integra- 
tions were carried out using the hybrid scheme w hich is 
part of the MERCURY package (jChambersI [i~999l) . This 
scheme is based on a second-order mixed variable symplec- 
tic (MVS) algorithm; it switches to a Bulirsch-Stoer scheme 
within a certain distance from a massive object. For all the 
integrations reported here, this distance was 2 Hill radii. 
The integration time step was chosen to be 4 days (l/20th 
of the orbital period of Mercury). During the integration, 
MERCURY detected and recorded close approaches with 
the terrestrial planets. We caution that, despite the time 
reversibility of the equations of motion in the absence of 



http:/ /hamilton. dm. unipi.it/astdys/ 
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Fig. 5. Minimum distances during encounters with Mars (left), Earth (middle) and Venus (right) as recorded by 
MERCURY during the backward (top) and forward (bottom) integrations of 1000 clones of (190491) 2000 FJ 10 . 



dissipation, backwards integrations will not, in general, pro- 
vide information of the dynamical history of the asteroid; 
instead, the forward evolution sta rting from the pos sible 
source regions must be considered (jBottke et al.ll2002l ). As 
we show in the following, however, the case of (190491) 
2000 FJio can be considered e xceptional in this sense. In 
addition, the Yarkovsky effect (Bottke et al. 2006) was not 
taken into account; we expect that its contribution to the 
dynamical evolution of the asteroid will be negligible com- 
pared to those of planetary close encounters and secular 
perturbations. 

5.2. Dynamical Evolution 

Figures [5j M an d Table [5] show the distribution of the close 
approach distances of the clones to Mars, Earth and Venus 
as a function of time, the distributions of their perihelia 
and aphelia and their orbital classification statistics. Note 
that we have used the Minor Planet Center definitions 
for the Amor, Apollo and Aten NEA classifications. These 
are slightly differe nt than those used by other authors (eg 
iBottke et al]l2002ft . 

At present the clones of this Amor asteroid are expe- 
riencing close approaches with Mars but not with either 
Earth or Venus (left panels of Fig. [SJ. Despite the clones' 
nominal perihelion distance being ~1.01 AU, this situation 
persists ±10 kyr into the integrations. At that time, es- 
sentially all clones have q > 1 AU (left panels of Fig. [5]) 
and remain classified as Amors (Table [5]). One likely con- 
tributing factor is a 3:2 near- resonance between the orbital 
periods of 190491 and the Earth, which renders close en- 
counters with the planet infrequent if not impossible (see 
Section E}. 

In the backward integrations, close approaches of clones 
to the Earth occur at a gradually increasing rate (upper 
middle panel of Fig. [5]) during a period of ~300 kyr. Their 



perihelia remain at or above 1 AU at —100 kyr (upper mid- 
dle panel of Fig. [5]); consequently they retain their Amor 
classification (Table [5]). By that time, only 15 clones have 
experienced close encounters with the Earth. In the forward 
integrations, the onset of Earth encounters is abrupt and 
occurs between +50 and +100 kyr. At +100 kyr most of 
the clones have perihelia < 1 AU and have become Apollos 
(lower middle panel of Fig. [5] Table [5]). This difference be- 
tween the outcomes of the forward and backward integra- 
tions is likely due to the secular evolution of the eccentricity. 
Fig. [7] shows the time series of a and e ±2 x 10 5 yr from 
the present for one in every 20 clones. Although the evo- 
lution of individual clones is inherently chaotic, we note a 
statistical trend towards smaller values of e in the past and 
larger values of e in the future. The timescale of variation 
(50+ kyr) leads us to suspect that it is related to one or 
more of the secular eigenmodes of the solar system ()Laskarl 
Il990f ) . The "critical" value of e required to reduce q below 1 
AU is 0.25 for a mean a value of 1.33 AU. Once e increases 
past that critical value - indicated by the gray dashed line - 
in the forward integrations, Earth encounters become pos- 
sible (see Fig. [5]) and, as a result, the scatter in a and e 
(bottom panels) increases significantly. In the backward in- 
tegrations, the critical value for e is never reached for all 
but a few cases. 

This result likely reflects a predictable event of the aster- 
oid's dynamical evolution. In other words the real asteroid 
will begin encountering the Earth sometime in the interval 
50-100 ky in the future and, as a consequence, become an 
Apollo asteroid. 

The final orbital distribution of the clones is shown in 
the right panels of Fig. [6] and the last column of Table [2 
Assuming Poissonian statistics (ie a = VN), the asteroid 
is more likely to be an Amor than not an Amor (3-tr level) 
at t = — 1 Myr. In the forward integrations, we can make 
the weaker statement that the most likely out of the four 
possible outcomes we examined is an Apollo classification. 
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Fig. 6. Distribution of clone perihelion and aphelion distances of 1000 clones of (190491) 2000 FJio at three different 
times during the backward (top) and forward (bottom) integrations. 



However, this result is only significant at the 2-a level when 
compared to that of an Amor classification. The final per- 
ihelia and aphelia of the clones appear to be more widely 
dispersed for the forward than for the backward integra- 
tions. Finally, we note that the forward integrations pro- 
duced a significant fraction of Atens and that two of the 
clones collided with t he Sun, a common en d result of NEA 
dynamical evolution (jFarinella et al.lfl994|l . 

5.3. Probability of collision with the Earth 

The distribution of close approaches as a function of dis- 
tance can be used to estimate the collision probability with 
the Earth. To do this, we first converted the recorded dis- 
tances qi and speeds v.; at closest approach to impact pa- 
rameters bi and velocities at infinity Voo,i through the ex- 
pressions (|Oniklll97rlD : 

Voo,i = \J vf - (1) 

bi = qJl + -^- (2) 

where [i denotes the product of the planet's mass with the 
gravitational constant G. 

Then we collected the bi into bins of width Abi = 0.005 
in units of the Earth's Hill radius (~ 0.01 AU) and fit- 
ted linear laws to the data for the backward and forward 
integrations separately, expecting that 

N&) = N(h < b < h + Ab>) ~ hAh. (3) 




itc -i 

0.00 0.05 0.10 0.15 0.20 0.25 



Impact parameter (Hill radii) 

Fig. 8. Linear fit to the number of close encounters to the 
Earth as a function of the impact parameter b for the for- 
ward (upper curve) and the backward (lower curve) inte- 
grations. 



Figure [8] shows that the frequency of deep close encounters 
(<C IRh ) is higher in the future than in the past, consis- 
tent with the higher dispersion of future over past clones as 
found earlier in this Section. The expected number of col- 
lisions with the Earth per Myr is N c = f Q E N(b)db where 
&e is calculated from Eq. [2] by setting q = Re, Earth's 
radius and assuming Voo^ =< Voo^ >~ 5.8 km s . This 
evaluates to ~ 0.11 and ~ 0.08 for the forward and back- 
ward integrations respectively. Hence the probability of im- 
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Fig. 7. Orbital evolution of asteroid clones 2 x 10 5 yr backwards (top) and forwards (bottom) from the present. The left 
panels show the semimajor axis while the right panels show the eccentricity. Only one in every twenty clones is plotted 
to maintain clarity. The grey dashed line denotes the minimum eccentricity required for an orbit with a semimajor axis 
of 1.33 AU to cross that of the Earth. 



pact of (190491) 2000 FJi with the Earth in the interval 
[-10 6 ,10 6 ] yr is - 2 x 10~ 4 . 

6. Synthesis: The Origin of 2000 FJi 

iBottke et all (|2002ft found that - 85% of all NEOs with 
H < 22 originate in either the inner (1MB; a < 2.5 AU) 
or central (CMB; 2.5 AU < a < 2.8 AU) Main Belt. 
About half (~ 53%) of existing Amors originate in the in- 
ner Belt alone. They are delivered in the terrestrial planet 
region via the 3:1 mean motion resonance with Jupiter, the 
v§ s ecular resonance a nd close encounters with Mars (see 
also lBinzel et al.| [2^004V The S-complex taxonomic classifi- 
cation deduced from our observations is consistent with this 
premise. The dynamical simulations show that it has been 
an Amor for at least 100 kyr from the present. The gradual 
loss of determinacy in the eccentricity evolution (Fig. [7]) al- 
lows us to extend the validity of this statement for up to a 
few 100s of kyr in the past. Since only a few clones attained 
perihelia > 1.3 au at the end of both the forward and back- 
ward integrations, we conclude that, if it arrived from the 
asteroid belt, it likely did so > 1 Myr ago. A small but sta- 
tistically significant difference between the outcomes of the 
full 1 Myr runs in the past and in the future suggests that 
the real asteroid is currently evolving from the Amor to 
the Apollo dynamical class. However, given the caveats in 
interpreting backwards integrations mentioned in Section 
5.1, we believe this conclusion to be tentative. To better 
quantify the likelihood of different scenarios, it would be 



necess ary to apply methods such as that of IBottke et al.l 
(l200l . 

We also note that > 95% of clones in both forward and 

o 

backward integrations maintained an inclination of < 15 , 
indicating that it is unlikely to have originated in the high 
inclination Hungaria and Phocaea families. If it originated 
within either the 1MB or the CMB (sour ce regions of 75% 
of Amors according to IBottke et al.ll2002[ ) it may be a for- 
mer member of a family dominated by S-type asteroids, the 
most populous of wh ich are the Flora and Eunomia families 
(|Zappala et al.lll995l) . 



7. Accessibility from the Earth 

To quantify the accessibility of the asteroid from the Earth, 
we constructed two way (Earth - Asteroid - Ear th) Type II 
kepler ian arcs based on the Gauss method as in iBate et ail 
(|197lD . Arrival and departure dates were determined by 
minimising the AT^ at the asteroid. We considered that a 
launch/return window existed when this quantity was equal 
to or less than that expended by the NEAR spacecraft in 
1999 to rendezvous with 433 Eros in early 2000 (-0.965 
km s -1 : iDunham et aT1l2002[ ). We found that consecutive 
launch windows for the asteroid are spaced ~3 years apart 
(upper panel of Fig. [5]) and that the same holds for return 
windows (lower panel of Fig. [9]) . Typical one-way trip times 
are ~1 year with a one-year wait at the asteroid before 
insertion into the Earth return trajectory. 
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Fig. 10. Round trip to 190491 using Type II two-impulse transfers. Left: outbound leg. The locations of the asteroid 
at departure and the Earth at arrival are denoted by the red disks. Right: Return leg. The red disks now indicate the 
location of the Earth at departure and that of the asteroid at Earth arrival. 
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speeds with 190491 at arrival. The AV expended by NEAR 
upon arrival at 433 Eros is indicated. Bottom: Return leg. 



to note for the discussion that follows that departure from 
and return to the Earth occurs when the asteroid is nearby. 

Interestingly, during the period 2000-2100 AD 
favourable launch windows occur in two groups, the first 
in 2000-2012 and the second in 2047-2059. Corresponding 
favourable return windows span the periods 2010-2020 
and 2058-2070 respectively. On those occasions rendezvous 
with the asteroid requires < 1 km s _1 of Ay and a similar 
amount for departure and return to the Earth. The total 
AV for arrival at, and departure from, the asteroid is < 2 
km s _1 on three round-trip opportunities in the period 
2052-2061. 

We find that this is the result of a 3:2 near-resonance 
between the orbital period of the asteroid and that of the 
Earth. This is best demonstrated if one views the aster- 
oid's trajectory in a frame that rotates with the Earth's 
mean angular velocity around the Sun (Fig. Ill[) . Due to 
the near-resonance the asteroid traces out a pattern with re- 
spect to the Earth completing one revolution every 3 Earth 
years (left panel) . The two loops correspond to the asteroid 
passing through the pericentre of its orbit. As of 2011 (left 
panel), one of those loops lies close to the Earth. Because 
the resonance is not exact, the pattern slowly precesses 
in a clockwise direction so that, after 15 years (i.e. 2026) 
close approaches to the Earth are no longer possible (mid- 
dle panel). Half a precession period later (18,000 d or 50 yr) 
the original configuration is recovered and close approaches 
to the Earth become possible again (right panel). In fact, 
querying the asteroid's ephemeris using the MPC online 
too{3 shows that the two closest approaches of the asteroid 
to the Earth (A ~ 0.1 AU) for the remainder of the 21st 
century occur in 2058 and 2061. 



A set of example trajectories generated by our code are 
shown in Fig. 1101 Departure from Earth occurs in Autumn 
2008 (voo = 5.5 km s _1 ) and arrival the following Autumn 
(AV = 0.95 km s _1 ). Departure one year later (AV = 
0.75 km s _1 ) brings the hypothetical spacecraft back to the 
Earth in October 2011 (v^ = 6.3 km s _1 ). It is important 



8. Conclusions 

The Amor NEA (190491) 2000 FJi is one of the most ac- 
cessible spacecraft targets. Its effective diameter of D e g — 
0.13±0.02 km places it in the transition zone between grav- 

3 http:/ /www. minorplanetcenter.net/iau/MPEph/MPEph. html 
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Fig. 11. The trajectory of the asteroid in a frame rotating with the Earth's mean motion over periods of 3yr (left panel), 
15 yr (middle panel) and 50 yr (right panel). The small ellipse at (1,0) represents the trajectory of the Earth in this 
frame, due to its orbital eccentricity (~ 0.017). 



itationally bound rubble-piles and monolithic bodies, he ld 
together by their internal strength ()Asphaug et al.|[2002f ). 

2000 FJio belongs to the S-complex of evolved asteroids, 
making it scientifically less interesting than the primitive 
objects that are usually the targets of robotic sample-return 
missions. On the other hand, its relatively large size, slow 
rotation, and accessibility from the Earth can make it a 
source of minerals and elements import ant for industry a nd 
a suitable target for a piloted mission (jAbell et al.l l2009). 

Due to its Earth-like orbit, during the next 10 years 
2000 FJio will make several close approaches to the Earth 
which can be used to determine its shape and spin axis 
orientation. Until 2020, there will be favourable apparitions 
of decreasing quality every 3 years when the asteroid will 
be brighter than V = 21 mag. This is rather uncommon 
among NEAs. 
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